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■ Abstract 
> , 

I A novel method for the exact solvabihty of quantum systems is dis- 

CN ' cussed and used to obtain closed analytical expressions in arbitrary di- 
mensions for the exact solutions of the hydrogenic atom in the external 

I potential AV{r) = 6r + cr^, which is based on the recently introduced 

' supersymmetric perturbation theory. 

I One of the challenging problems in nonrelativistic quantum mechanics is to 

find exact solutions to the Schrodinger equation for potentials that are used in 
^ ■ different fields of physics. In particular, the perturbed Coulomb potentials repre- 

sent simplified models of many situations found in atomic, molecular, condensed 
matter and particle physics. There has been much interest in obtaining analytical 
solutions of such potentials in arbitrary dimensions. These problems have been 
studied for years and a general solution has not yet been found. 



' Such class of potentials, 

V{r) = —a/r + br + cr'^ , (1) 

are possible candidates for the quarkonium potential as has been indicated by the 
quarkonium spectroscopy. In the special case of c = and 6 > such potentials 
reduce to the well known charmonium potential. Apart from its relevance in 
heavy quarkonium spectroscopy, this class of potentials with c = has important 
applications in atomic physics. The Stark effect in a hydrogen atom in one 
dimension is given exactly by the charmonium like potential (6 being the electric 
field parameter). The more general class of these potentials with c > is also 
relevant in atomic physics. This could be interpreted as the potential seen by 
an electron of an atom exposed to a suitable admixture of electric and magnetic 
fields. In addition, nuclei in the presence of an electron background form a 
system which is important for condensed matter physics and for laboratory and 
stellar plasmas. The potential between two nuclei embedded in such a plasma is 
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approximately Coulomb plus harmonic oscillator, which corresponds to 6 = in 



As the exact form of such interactions are being unknown to a great extent, 
it is thus desirable to study the general analytical properties of a large class of 
potentials in (1). In connection with this, the analyticity of the energy levels 
for these kind of potentials was investigated rigorously by many authors using 
different theories [1]-[14] in relation to their potential applications in spectroscopic 
problems. 

In this letter, we introduce an alternative, simple formalism for an algebraic 
solution of the Schrodinger equation with the perturbed Coulomb potential and 
find exact solutions in A^-dimesional space. The new formalism is based on the 
supersymmetric quantum mechanics and uses the spirit of perturbation theory, 
in which there is no room for the drawbacks encountered in the calculations with 
the perturbation theories available in the literature. The work introduced in this 
letter also clarifies, within the powerful frame of the present formalism, that the 
potential in (1) is in fact an exactly solvable shape invariant potential unlike 
the claim in a recent work [11] where the authors stated that the supersymmetric 
quantum mechanics yields exact solutions for a single state only for the perturbed 
Coulomb potential in (1) which was treated as a quasi-exactly solvable potential. 

Let us now introduce the formalism. The reduced radial wave Schrodinger 
equation for a spherically symmetric potential in A'"-dimensional space reads 



where Vq is one of the exactly solvable potential, such as the Coulomb and har- 
monic oscillator potentials used through the present work, n = 0, 1, 2, . . . being 
the radial quantum number, and A = (M — 3)/2 with M = N + 2£. We see 
that the radial Schrodinger equation in N dimensions has the same form as the 
three-dimensional one. Consequently, given that the potential has the same form 
in any dimension, the solution in three dimensions can be used to obtain the 
solution in arbitrary dimensions simply by the use of the substitution i ^ A. In 
the above equation, AV stands for the small perturbation (= br + cr^) assuming 
the potential parameter a is large when compared to other coupling parameters 
b and c in (1) with the consideration of perturbed Coulomb interactions while for 
the case of large c and relatively small a and b the perturbing potential involves 
AV = ~a/r + br leading to a perturbed harmonic oscillator potential. This 
significant feature will be used later in this work to test our results. 

Using the spirit of the standard perturbation theory, we put forward here 
a more general formalism to treat quantum mechanical perturbation problems 
efficiently in a simple framework, which has been already discussed in detail 
through the recent works [15, 16] involving different applications of the model. 
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y(r)-K , V{r) = VEs{r)+AV{r) , VEs{r) = Vo{r)^ 



A (A + 1) 



2mr^ 



(2) 
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Start with 



*n('^) = Xn{r)(l)n{r) , 



(3) 



in which Xn is the known normahzed eigenfunction for the exactly solvable po- 
tential in (1) including the barrier term and <;/)„ is a moderating function due to 
the perturbing potential. Substituting (3) into (2) yields 



n Xn 



(4) 



With the new definitions, 



one arrives at 



2m Xn 



^ Xn 



'2m Xn 
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(5) 
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where is the eigenvalue for the exactly solvable potential of interest, and 



2m I 0n Xn 4>n, 



/^Wl - -j=^W'^ + 2W^/^W^ = A1/ - Ae„ , (7) 
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in which Ae„ is the correction term to the energy due to Ay, and — e„-|- Ae„. 
Subsequently, Eq. (4) reduces to 



h 



'2m 



{Wn + AWn)' = V-En • 



(8) 



As the complete spectrum and wave functions are known in the analj^ical 
form for the exactly solvable potentials appeared in (6), one needs here to solve 
(7) to obtain in a closed form the corrections to both, energy and wave function. 

Proceeding with the perturbed Coulomb potential in arbitrary dimensions. 



V{r) 



_a ^ A(A+ 1) 
r 2mr^ 



+ br + cr^. 



(9) 



and set the superpotential, for the ground solutions, 

iA + i)h 



'2mr 



(10) 



leading to the first part of the potential in the bracket, and from the literature, 
corresponding normalized wave function and energy are 



Xn=o{r) = iVor'^+^exp 



ma 



(A + 1)^2 
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ma 



en=o 



2r (A + 1)^ 



For the perturbing potential, AV — br + cr^, the physically acceptable unique 
choice is 

AWn=o{r) = V~cr , (12) 
which satisfies Eq. (7) from where one readily sees that 



Ae„=o = , & = 



Ama 



{M-l)h 



(13) 



It is stressed that such solutions in general have constraints on the potential 
parameters as appeared in (13). These constraints differ for each eigenvalue, 
and hence various solutions do not correspond to the same potential and are not 
orthogonal. 

Prom (12), the moderating function 



0n=o(r) = exp 



2m 



h 



J AWn=o{r) 



— exp 



b{M-l) 



(14) 



Therefore, the full wave function ^n=o for the potential in (9) takes the form 

2ma 



*„=o(r) = Xn(r)(/)„=o(r) ^ No 2 ) exp 
Finally, the exact ground state energy is 

En=0 — ^n=0 + Aen=0 - 



(M - 1) h- 



6(M-1) 
' 4a 



2ma 



r {M-iy 



+ 



M (M - 1) bh^ 
Ama 



(15) 
(16) 



The results in Eqs. (15) and (16) are exact and agree with those in [11, 12]. 

As noted earlier, the potential in (9) behaves also like a harmonic oscillator in 
case for large c and relatively small a and b. This enables one to check out explic- 
itly the results obtained above and the reliability of the formalism introcuded. 
For this consideration, we set the superpotentials 



v2mr 



(17) 



which yield, through the use of either Eqs. (6) and (7) together, or (8) alone, 

e„=o = (2A + 3) , Ae„=o = , (18) 



which are exactly equivalent to those found for the perturbed Coulomb case in 
(16). Similarly, using (17) one readily arrives at (15). 
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For the generalization, we further make clear that the present technique is 
also applicable to excited states, for which one needs to use the shape invariance 
property and the relation between supersymmetric partners [17] 

V+ir, ao) = V'ir, ai) + R{ai) , (19) 

where ao = A and ai = f (ao) = A + 1 are the position independent parameters 
while V"^ are the supersymmetric partners 

{Wn=o + ^Wn=of T {Wn=o + ATK=o)' ^ - E'^^ , (20) 

v2m 

which depend upon superpotentials in (10) and (12) for the case a > b,c. Bearing 
in mind that the replacement of A with A + 1 in V~{r,A) leads to V~{r,A + 1) 
and keeping the potential parameters b and c fixed as in [11, 12] while allowing 
the change in a, see Eq. 13, we obtain by the straightforward calculations 

= + + 

K = i:it:(a.) = -^[(2A + 3)-2(n + A) + 3] , 
k=i v2m 

K = _^ + ^[2(n + A) + 3] , 71 = 0,1,2,... , (21) 



which agrees with [12] in A^-dimensional space and with [9] in three dimensions. 
Furthermore, one can easily construct the bound n*'* state wave functions from 
the ground state wave function using the supersymmetric definition [17] 



^„+i(r,ao) oc A+(r,ao)^„(r,ai) , A+(r,ao) = 



^ <=o(r,ao) 



dr^ *n=o(?^, ao) 



^/2m 

(22) 

We conclude with two remarks. First, the present formalism can be general- 
ized to all the polynomial forces 

Kn(r) = Ar^"" + Br^""-' + ... + F/r + G/r^ (23) 

as an alternative but simple treatment to the other works [18] , and the references 
therein. With the consideration of the perturbed Coulomb problem we have 
clarifed that the explicit solution of the related Schrodingcr equation with the 
potential family in (23) remains feasible in an almost complete parallel with their 
m — 1 predecessors. Along this line, the works are in progress. 

The second remark we wish to make is that apart from the inherent interest 
one has in the existence of exact solutions, the results reported here are likely 
to be useful in perturbation calculations for the excited state energies and wave 
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functions, particularly if the method of other perturbation theories appeared in 
the literature for evaluating second- and higher-order corrections can be extended 
to cover such cases. In this context, we believe that the simple form of our new 
approach to such problems opens a new direction of development towards many 
practical applications yet to be constructed and appreciated. 
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